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Abstract. We discuss the various contributions to the jet luminosity 
in proton blazar models of active galactic nuclei and describe a method 
of estimating the jet luminosity from the observed spectral energy dis- 
tribution (SED) and the fitted model parameters. We apply this to a 
synchrotron proton blazar (SPB) model for Markarian 501. 



1. Introduction 

The mechanical lumino sity of an AGN jet of cross-sectional area A with Lorentz 
factor r = I/a/I — /3 2 , containing jet-frame (primed variables) matter density 
density p', energy density u' and pressure p' is given by 

£jet,mech = T 2 p C A[p'c 2 (T - 1)/T + u' + p'} . 

(see e.g. Leahy 1991). To calculate the total jet luminosity Lj et , measured in the 
rest frame of the galaxy, we adapt the formulae of Bicknell (1994) and Bicknell 
&; Dopita (1997), given for the synchrotron self-Compton model, to apply for 
the case of proton blazar models. In this paper we shall apply the formula 
adapted for proton blazar models to a synchrotron proton blazar (SPB) model 
for Markarian 501 described by Miicke & Protheroe (2000). 



2. Proton blazar models 

In proton blazar models, the high energy part of the SED is due to interaction 
of protons accelerated along with electrons in the AGN jet. The interactions are 
pion photoproduction on either low energy photons of the low energy part of the 
SED produced as synchrotron radiation by electrons in the jet (e.g. Mannheim 
1993), or on direct or scattered UV bump radiation from an accretion disk (e.g. 
Protheroe 1997), and direct synchrotron emission by protons, muons and charged 
pions (e.g. Miicke & Protheroe 2000). To accelerate protons to sufficiently high 
energies that they can produce the high energy part of the SED, a large magnetic 
field is required. 

Proton blazars would contain relativistic plasma of electrons and protons, 
and a non-negligible magnetic field. For a n' e (^ e ') = n e0 ^ e '~ 2 power-law (7^ < 
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7e' < le'2) f° r the jet-frame number density of electrons, the total number density 
of (relativistic) electrons is n' e = w e o/7ei) and the energy density of (relativistic) 
electrons is 3p' e = u' e = n' e m e c 2 {y e ') = n e0 m e c 2 ln(7 e / 2 /7 e / 1 ). Similarly, for a 
n 'p(lp') = n polp'~ 2 power-law (7^ < 7 p ' < 7 P2 ) for the jet-frame number density 
of protons, the total number density of (relativistic) protons is n' p = « P q/7 p1 , 
and the energy density of (relativistic) protons is 3p' p = u' p = n' p m p c 2 (7/} = 
n P0 m p c 2 In(7p' 2 /7 P ' 1 ), giving 

/ 2 _ 3 Pe / 2 _ 3 ^P 
n m e c — — - — — — p-, n m p c — — - — — - — — -. 

7e'l ln(7e 2 /7ei) p 7/1 Hl P2 nPi) 

Assuming that the number of relativistic electrons will be greater than the num- 
ber of relativistic protons, and applying charge conservation, i.e. the number of 
'cold' (non-relativistic) protons equals the number of 'hot' (relativistic) electrons 
minus the number of hot protons, one obtains 



-^jet, mech — 

= T 2 l3cA 



m r . 



(T-l) 



m e 7 e i ln(7 e / 2 /7 e / 1 ) j p [ ln(7 p / 2 /7 p / 1 ) 

¥ P r 2 (3cA [x P ^^ + 1 + P J + 4 

Pv 



+ Ap' e + Ap' p + Ap' B 



Pp fp 

since for a tangled magnetic field u' B = 3p' B , and where 
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m e p' p 7^ ln(7 e2 /7 e / 1 ) 7^ ln(7p / 2 /7 P ' 1 ) 



We now consider how to apply these formulae to proton blazar models. First, 
consider the radiation efficiency of protons which depends on their effective 
energy loss rates r' for synchrotron radiation (syn), photoproduction leading to 
electromagnetic radiation (EM) or anything (any), and adiabatic losses (adiab), 
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Averaging over the input energy spectrum gives the total radiation efficiency 



/ lp'lp~ 2 djp' 



For electrons, we assume the radiation efficiency to be ( e = 1. 

If we assume that the low energy part of the observed SED is due to elec- 
trons, and the high energy part is due to protons, then we may infer values of p' e 
and p'p directly from the observed bolometric flux in the two parts of the SED 
for an assumed Doppler parameter and emission region size. The bolometric 
luminosities for the two parts of the SED (in any frame since they are Lorentz 
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invariant) are then L^J = ( e r 2 (3cA4p' e and LjjjJ 1 = C p T 2 (3cAAp' p . These bolo- 



metric luminosities are related to the observed bolometric fluxes (corrected for 
relativistic beaming) from the two parts of the SED by L l g£ = And 2 - S l ° w 1 n2 



and L^f 

d,L the source's luminosity distance. Hence, we obtain 
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4Trd 2 L S^/D 2 with D = [T(l - PcosO)}- 1 the Doppler factor and 
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where p' B = [B' 2 /(2/x )]/3, and 
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3. Application to SPB model for Markarian 501 

We use these formula to estimate the total jet power of Mrk 501 during its 1997 
flare to be ~ 10 46 erg/ s, and find the contributions to the total jet power of cold 
protons, magnetic field, and accelerated electrons, relative to that of accelerated 
protons. Fig. 1 shows the dependence of the total jet luminosity on the Doppler 
factor D for a fixed variability time scale i var = 12 hours, and with B, n p ' and 
a target photon density appropriate to Mrk 501 during flaring. Clearly visible 
is the fact that in hadronic models the proton kinetic energy and the poynting 
flux dominate the total jet luminosity, while the electron kinetic energy is only 
of minor importance. At high Doppler factors the emission region becomes so 
large that one needs only relatively small magnetic fields and proton densities 
to fit the observations. In addition, adiabatic losses become small resulting in a 
decrease of the required kinetic proton luminosity. For example, B ~ 5 G and 
n' p ~ 10~ 2 cm~ 3 are sufficient to fit the Mrk 501 flare for D = 50, while for 
D = 8 magnetic fields of over 30 G and proton densities of n' p w 10 4 cm~ 3 are 
needed. The total jet luminosity exhibits a minimum of ~ 10 46 erg/s atD« 12. 

In the framework of the jet-disk symbiosis (e.g. Falcke & Biermann 1995), 
the jet luminosity should not exceed the total accretion power Q& CCT for the 
equilibrium state. Accretion theory relates the disk luminosity to the accretion 
power. Page & Thorne (1974) give -Ldisk ~ (0.05 — 0.3)Qaccr- Disk luminosities 
for 'typical' radio-loud AGN lie in the range Ldisk ~ 10 44 — 10 48 erg/s with BL Lac 
objects tending to the lower end on average. Specifically, for Mrk 501 there are 
no emission line measurements available, and this complicates the evaluation 
of its disk luminosity. However, any observed UV-emission in the flaring stage 
may put an upper limit on it. Historical data give -Ldisk ~ 10 43 — 10 44 erg/s 
(Mufson et al 1984, Pian et al 1998), and we obtain for the accretion power, 
Qaccr ~ (3 — 200) x 10 43 erg/s, at least a factor 5 below the value necessary 
to comply with the constraint of the disk-jet symbiosis. Note, however, that 
the estimate of Qaccr is based on archival non-flaring data from Mrk 501, and 
we could speculate that either the disk has pushed more energy into the jet 
during TeV-flaring, or that the flaring stage can not be considered as a steady 
state. Also, accretion theory might predict larger conversion efficiencies of the 
accretion power into disk radiation than actually might occur in BL Lac objects. 
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Figure 1. Dependence of total jet luminosity Lj et and its contribu- 
tions (Lb - magnetic field, £ Pi hot an d -^ P , C oid ~ hot and cold protons, 
L e - electrons) on D for model parameters which reasonably fit the 
Mrk 501 1997 flare SED (see Miicke & Protheroe 2000 for details). 
L 40 = L/10 40 erg/s. 
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